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Abstract 

An Enskog-Landau kinetic equation for a many-component system of 
charged hard spheres is proposed. It has been obtained from the Liouville 
equation with modified boundary conditions by the method of nonequihb- 
rium statistical operator. On the basis of this equation the normal solutions 
and transport coefficients such as bulk k, and shear 77 viscosities, thermal 
conductivity A, mutual diffusion D"^ and thermal diffusion have been 
obtained for a binary mixture in the first approximation using the Chapman- 
Enskog method. Numerical calculations of k, 77, D"^ and for mixtures 
Ar-Kr, Ar-Xe, Kr-Xe with different concentrations of compounds have been 
evaluated for the cases of absence and presence of long-range Coulomb inter- 
actions. The results are compared with those obtained from other theories 
and experiment. 

PACS: 05.60.-Hw, 05.70.Ln, 05.20.Dd, 52.25.Dg, 52.25.Fi. 
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1. Introduction 

The Enskog-Landau kinetic equation for an one-component system of charged hard 
spheres |l) is comparatively new. Being obtained from the first principles of statis- 
tical mechanics it describes well transport processes in moderately dense systems. 
Its collision integral does not have any singularities on short-range distances. Nu- 
merical calculations of shear viscosity rj and thermal conductivity A on the basis 
of normal solutions to this equation using the Chapman-Enskog method are in a 
good agreement between theory and experiment . In next papers ^ , nonsta- 
tionary solutions and transport coefficients have been obtained using the boundary 
conditions method. In the limiting case of stationary processes, the results of |2j 
and iQ] coincide between themselves. Despite this fact, the one-component system 
remains always by a model system [f[ 0, §■ 
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Traditional interest in many-component systems, including two-component mix- 
tures of dense gases, liquids and plasmas, has been intensified recently by the con- 
siderable advancement in computer equipment. This allows to supplement theory 
by results of mathematical simulations rather easily. 

The transport coefficients of nonequilibrium systems have been investigated us- 



ing different methods, namely: Bhatnagar-Gross-Kruk theory Green-Kubo |fO 



and Frost methods, at unified description of kinetics and equations of non- 
linear hydrodynamics [^ , in one-liquid approximation Mathematical con- 
struction of extended thermodynamics of dense gases and liquids (including hard 
sphere mixtures within RET theory [|l|) has been carried out jl4|, thermodynami- 
cal compatibility of Enskog-type kinetic equations for M-component systems and 
Enskog-like kinetic equations for reacting mixtures has been proved |p^ , results 
of macroscopic and microscopic theories have been compared [0. At the same 
time, when methods for calculating transport coefficients were improved, numer- 
ical simulations such as Monte Carlo (MC) and molecular dynamics (MD) were 
optimized |l^ also. In some cases, purely theoretical problems lead to some 
other ones: the Gross- Jackson model [pO[, inverse problem of thermal conductivity 



[ |2l[ , effective transport coefficients 22, |2^, In such a way the efforts to ex- 
pand the existing theory || on many-component systems looks quite naturally. 
Moreover, a sequential theory for charged particles does not exist yet. In papers 
|, [10[ |ll|, |l2[ 1^, 1^, |2l], m 0, transport coefficients 



(such as viscosity rj, conductivity cr, thermal conductivity A, mutual diffusion D"^ 
and thermal diffusion D^) have been calculated in one or other way. It is pointed 
everywhere, that the best coincidence of theoretical and experimental data (as well 
as MC or MD) is obtained for neutral systems at comparatively not high densities. 
As was noted in the coincidence at high densities is bad (theory does not 
work). It is summarized in on the basis of MC simulations, that rj coincides 
with experimental value better than A in major cases. Generally, a worse agree- 
ment of theory and experiment for A is known for a long time. First of all, it is 
connected with a strong growing of interparticle correlations as density decreases. 
These correlations should be taken into account more precisely . In its turn, this 
leads to changes in local conservation laws, which are used at finding for solutions 
to kinetic equations 

In the present paper, the Enskog-Landau kinetic equation for a multicomponent 
system of charged hard spheres, its normal solution by means of the Chapman- 
Enskog method, and the calculation of transport coefficients for binary mixtures 
are proposed. The theoretical part of this work was published for the first time 
in papers [p6| |2^ . Now it is upgraded by some results of numerical calculations 
for transport coefficients as well as their comparison with some experiments and 
MD data. A special attention is for i7-integrals calculation for long-range type of 
interaction (Coulomb interaction in our case) and for finding of static screening 
radius. Expressions for spatially dependent distribution functions in relations for 
transport coefficients are discussed concisely as well. 



2 



2. Initial relations 



Let us consider a nonequilibrium system of N classical particles, which consists of 
M kinds of charged hard spheres. Each kind has Na particles, where a is varied 
between 1 and M, with X]a=i — N. The masses of particles are equal to 
uia, altogether they occupy a volume V. It is agreed that the thermodynamic 
limit transition N oo, V — > oo, N/V — const takes place. The charges of 
particles are equal to ZaC^ where e is the electron charge, Z^, G Z. The condition 
of electroneutrality is applied to the system as a whole. It means that in the 
case of differently charged particles X]q=i ^a^a — 0, while in the opposite cases, 
where only positively charged particles are considered (mixtures of ionized gases, 
for example), a neutralizing continuum is introduced. 

The Hamiltonian of such a system can be presented as follows: 



a=lja = l a=l [i=l ja = lkf3 = l 

jai^kfj at a=/3 



where 

is the total interparticle interaction potential, which consists of the short-range 
(hard spheres) <P^^ and the long-range (Coulomb one in our case) interactions: 

<?'^^(ir-.„.,i) - lim ninM) - 1 o' l^i"! > 




A nonequilibrium state of this system is described by the nonequilibrium TV- 
particle distribution function g [x^; t) . Let us consider now much precisely the case 
of a two-component system, when M = 2. Then the nonequilibrium TV-particle 
distribution function, 

— Q {xi ,X2,---, Xjy^ , Xi, X2, ■ ■ ■ , Xjy^ ; t) , 

satisfies the Liouville equation and normalization condition 0, H. In the two- 
component case the iV-particle Liouville operator has a different structure, which 
takes into account a many-component nature of the system: 

a=lja = l a=ll3=ljo,=lkf, = l 

ja,9^kf3 a=f3 
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L{ja) = -i 



LUa, kp) = i 



dp. 



dPks 



3. Kinetic equation 



The kinetic equation for a-kind one-particle distribution function /i (a;" ; t) can be 
obtained similarly to an one-component system, but with some distinguishes, which 
are due to a multicomponent nature of the system under consideration. In such a 
way, in the "pair" collisions approximation this equation reads: 



' 13=1 

hm e'^^^g2{r'^,rl; t + r)h{x^;t + T)/i(a;f ; t + r), 



(3.1) 



where 



Equation is formally distinguished on analogous one for an one-component 

system by the presence of kind summations on its right-hand part in the colli- 
sion integral. The Enskog-Landau kinetic equation for a two-component system of 
charged hard spheres has been obtained for the first time in [M: 



d_ 

at 



ii(la) 



fiK;t) 



(3.2) 



4*'^ (x? ; t) + (x? ; t) + Imf (a:? ;t) + h (x? ; t) , 



where terms in the right-hand part are collision integrals. Explicit expressions for 
the collision integrals depend on a specific interparticle potential of interaction. 
The first and the second terms are collision integrals of the revised Enskog theory 
(RET) S|: 



M 

E 



I^"{x^;t) = / dt;f de dhhg'^P gf{cj^p\n,(3) x 

'/i(r?, ; t)h{rl ; t) - /i(r?, <; i)/i(r-^, v^^- 1) ] , 



(3.3) 



4'^K;i) 



M 

/3=1 



a0\ 
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(3.4) 



h{rl.<: t)Vh (r^ , ; t) - /i (r? , v", ; t) V/i (r^ , t;^ ; t) 



/3 



+ 



A (r?, t;f ; i)/i (rf , v^,'; t) + h (r?, i^?; t)/i (r^ , ; i) 



/3 



The next term is caused by the influence of long-range interactions in the mean 
field approximation (KMFT) |: 

(3.5) 



M 

— V 



1 



drf ^^'H|r?2''l)52^(r-?,rf;0^/i(x?;t)n^(r^;0, 



the last term is a Landau-like collision integral ^ : 



'9-' X 



(3.6) 



/3=1' 



(/i (r? , ; /i (rf , «f ; i) - /i « , ^? ; t)h (r-^ , ; 



This collision integral is presented here in a reduced Boltzmann-like form. It could 
be obtained going from the cartesian to cylindrical reference frame, entering im- 
pact parameter h, azimuthal angle of scattering e, distance along cylinder cen- 
treline ^ and integrating with respect to ^ and taking into account the condition 
g2{ri,r2\t) ^ 1. If one solves this kinetic equation using the Chapman-Enskog 
method, this form for the la test col lisio n integral is more convenient |Q. Other 
conventional designations in (|]^) ^ (3_^) are as follows: 



b 

e 

92 

m* 
rua 
n 

Ua 
-a/3 

e{x) 



impact parameter, 

inversed local temperature analogue, 

azimuthal angle of scattering, 

relative velocity of a- and /3-kind particles, 

two-particle correlation hmction, 

reduced mass, 

partial masses of particles, 

total density of particles number, 

partial densities of particles numbers, 

unit vector along r"^ direction, 

unit step function. 



An additional point for emphasizing is that vf and are hard spheres velocities 
after collision: 



12 



-Q/3/.a/3 c 

' 12 y 12 y 
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as well as v"* and V2* are velocities of charged particles after Coulomb scattering: 



(3.7) 



4. Normal solutions by means of the Chapman-Enskog 
method in the first approximation 

As usually, to a solve kinetic equation by means of the Chapman-Enskog method 
in the k th approximation we will use conservation laws for a set of hydrodynamical 
variables in the (/c — 1) th approximation. The components of the additive invariant 
vector could be chosen similarly to ^ , namely the hydrodynamical mass density, 
momentum and energy: 

M „ 

p{ri;t) = ^ dv^ fi{x'^;t)ma, 

a=l 
M » 

p(ri;t)y(ri;t) = ^ / d< t)m„<, 

a=l 

p(n;t)a;k(ri;t) = / d< A K ; t) ^^^^^ . 

The conservation law for p{ri;t) have the form of a continuity equation, whereas 
the equation of motion and equation of energy balance read as 

^^"■^'^^ dt ^ -a^9*(r-i;t)-P*(n;0:^V(ri;t) + 

where 

VJ^(ri; t) is the diffusion velocity of a-kind particles: 
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here (. . .) denotes an average value over the one-particle distribution function 
fi{xi;t). In the zeroth approximation, when the function /i(a;";f) is taken as 
a local-equilibrium Maxwell distribution function (k is the Boltzmann constant), 

, 2 ^ 



/f)(x?;t) = rvK;i) ( ) exp 



2nkT{rf;t)J "^1 2kT{rf;t) j' 
one obtains for P '(ri;t), q*(ri;t), V'^{ri;t) the following expressions: 

M 

pk^7pk^ pk^^n,A:T, (4.2) 



2 ^ 

P ' = 0, q^ = q^'' = = q' = 0, (^i! ^) «• 



Using the mass conservation law and equation (O) in the zeroth approximation 

and taking into consideration an explicit form for P '(ri;t), P *(ri;t), q*{ri;t) 
(4.2), we will find for the solution to equation ( [3.2| ) in the first approximation as 
the combination: 

/i(x?;t)-/f)(x?;i)(l + ^K;i)), (4.3) 

where the correction (p(x";t) is expressed through the Sonine-Laguerre polynomi- 
als: 



\2 



2 ' 



2fcT^ 2fcr M'^i'^)^'^^ 



2fcr \ 2fcr J ^ 3 

Here 

oo 

i?"(x) = ^i?:Lf/2(:r), 

oo oo 



n=0 



n!P(rt -t- r + 1) 



m\r{m -\- r -\- l)r{n — m + 1)' 
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d" = d"(r";t) is the diffusion thermodynamic force of a-kind of mixture 
Quantities satisfy the condition 



_:d"(r?;t) = 0. 

a = l 

In the case of an one-componen t sy stem, diffusion t herm odynamic forces are absent. 



If one puts the solution of ( [l.g| ) into equation (3^), one shall obtain a nonuni- 
form integral equation of Fredgolm-type. Therefore, additional conditions for the 
existence of solutions to this equation can be written down at once: this is or- 
thogonality of the right-hand side of (^^) up to solutions of the corresponding 
uniform equation. The last ones are parameters of the abbreviated hydrodynamic 
description. Whereas, additional conditions after mathematical transformations 
read: 

77i6 _ '^a r^a ™ Ab _ ""-a ma 

rib \ rrib Ub \ nib 

This should be taken into consideration any time one looks for expansion coefficients 
i?", A!^, B". In such a way, several independent infinite systems of equations 
appear. It is impossible, of course, to find an exact solution to these systems. One 
should restrict by a few first terms in expansion series. As a rule, practically, one 
considers quantities Eq, Bq and ^q, A" only. We should take in mind the fact, 
that the collision of identical particles gives no affect due to mass, momentum and 
energy conservation laws at elastic scattering. 

The basic relations to obtain A" and Af coefficients as well as their symbolic 



Appendix B 



solution are presented in Appendix A. The same for Bq coefficients is done in 



5. Transport coefficients 



Having one-particle distribution function (4.3) in the first approximation, we are 



able to write down the expression for mass density flow of a-component of mixture 

where D^f^ is the mutual diffusion coefficient, and is the thermal diffusion 
coefficient of a-kind of mixture. In the zeroth polynomial approximation 



mpn V 2ma 



kT 

D^ = man^\l- A^. (5.2) 

Then the expression for thermal diffusion ratio is easy to write as 
p Dt lA^ 



h 



T 
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Final expressions for coefficients of thermal diffusion one can write taking into 
account the structure of -coefficients (Appendix A.). To calculate quantities Eq 
it is convenient to consider new variables, namely, the variables of inertia centre. 
By that transformations these quantities are calculated exactly. This allows to 
write the expression for mutual diffusion coefficient in an analytical form: 



37rm/3 jTrm. 
Spric 



Stt / TrfcT 
8n 



Here "^f^l^^^^ "''j?!'''^-' are known as so-called f2-integrals |3^. The f2-integrals 
appear in calculation of collision integrals (integral parentheses [pO|) in kinetic 
theory of gases, liquids and plasmas. The explicit form of ]7-integrals depends on 
type of interparticle interaction, that is taken into consideration. In our case, the 
total interparticle interaction potential consists on hard sphere and Coulomb parts. 
In such a way, two types of i7-integrals will appear. 



a/3 q(i''P) 



ha 



dy e ^ y 



2p+3 a/3^(r) 



hs ' 



a/3 f^i-^.P) _ 



27r dbb 



alii 



X {b, y) 



a/3o('-) ^ 



2tt dbb 



1 — cos' 



Here x' is the scattering angle of hard spheres after collision, x* is the scattering 
angle of charged particles. The analytical calculation of /^-integrals presented in 
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Except the thermal and mutual diffusion coefficients in the first approximation, 
the stress tensor P '{ri;t) and heat flow vector q*{ri;t) one can calculate also and 
obtain other transport coefficients. The calculations give: 

P {ri;t) = P{ri;t) J - k{V : V{ri;t)) -2r]S {ri]t), 

M 

q{ri;t) = -AVT(ri;i) + ^ c^„rf". 

a=l 

Here k is the total bulk viscosity coefficient of mixture: 
8 m* 
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a,/3=l 



a/3, 



a, 13=1 



(5.3) 
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1] is the total shear viscosity coefficient: 

M 

5 ■ 2 



3 

77 = —K 



1 *^ 



a=l 
M 



15 

/3=1 

A is the total thermal conductivity coefficient 



A = - ^ fcm^ 



m/3 



a, 13=1 
M 



1 1 (nia - mpY 
-{rria + mp) - ■ 

z 8 nrio, + mp 



5 ^ , 2kT 



/3=1 



A" 



5 /2fc3r3 

4 V TO„ 



^0- 



(5.4) 



(5.5) 



There is important to consider transport coefficients (5.1) - ( |5.5| ) in a limiting 
transition, when the number of kinds M goes to unity. In such a case one must 
obtain results of theory for an one-component system of charged hard spheres in a 
neutralizing background j|] . This is easy to see if one puts 



ma ^ nib ^ m, 
Ha = nb = n/2, 



Bq - Bfj - Bo. 



By this means, there are no diffusion and thermal diffusion coefficients in an one- 
component system due to the absence of thermodynamic diffusion forces d". But 
the equality — appears as the Fredgolm condition. 



6. Numerical calculations 



In our paper numerical calculations have been carried out for transport coefficients 
Z)"'' (O), Dq, (O), K (U) and 77 (O). We have chosen mixtures of Ar-Kr, 



Ar-Xe, Kr-Xe as objects of our calculations. Neutral and charged mixtures of 
these gases are appropriate well for the model of neutral and charged hard spheres, 
which is considered here. This is so because other mixtures can have complex 
molecular structures. The transport coefficients have been calculated at several 
fixed dimensionless densities Z\ in a wide physically real for these mixtures ranges 
of temperature {A = X^qLi ^qO'q; = Trnc/B, cr is a hard sphere diameter). 

Magnitudes for neutral systems [Za = 0) of some transport coefficients {D"^ , 
rf) of Ar-Kr mixture have been compared with results of MD |32| and HSED 
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theory |3^] for two fixed temperature points Ti = 508 K and T2 — 1196 K at the 
same dimensionless density A = 0.1 {ua — np of course). The coincidence of these 
values is exact (see Table |^) . 

Table 1: The comparison of theoretical calculations of total shear viscosity and 
mutual diffusion D"^ with results of molecular dynamics simulations and other 
theories for Lenard-Jones liquids of argon and krypton. * - data was borrowed 
from [32], ** - data was borowed from [31], *** - data was borrowed from [33]. 



A 


T, K 


ri, Pa-s, MD 


■q, Pa-s, EDHST 


•q, Pa-s 


' S ' 


2^Ar — Kr 

' s 


0.1 


508 


6.47-10--"^ * 




2.58-10"^ 


1.518-10-^ 


2.568 10-'' 


0.1 


1196 


7.72-10"''^ " 


8.56-10-^ *" 


3.9510-^ 







Values of shear viscosity of argon component rjp^^ of neutral mixture have 
been compared with results of [^^ |3^, ^ . Theory parameters ctq, were taken 
fro m Ips], ^ Jiol, E2|. The obtained results are in a good agreement with 
[H, |35|, ^, |37|, |3^,]39[ |o[ |l], || data. All obtained temperature dependencies of 
all transport coefficients at the absence and presence of long-range Coulomb inter- 
particle interactions were published in our recent preprint [^ . But in the present 
paper we give the combination of some numbers of them only, which reflect good a 
general tendency. The difference up to two orders in magnitudes for the values of 
the presented transport coefficients is explained as follows: firstly, different kind of 
interactions; secondly, different densities. Both these factors affect on magnitudes. 
This testifies that the obtained formulae are very sensitive on change of system 
parameters. 

Graphs of transport coefhcients Z?"^, and r/ in the case of their temperature 
dependencies are shown in Figure 1, kind concentrations are the same [ua = n/s). 
Each graph is assembled by pairs: neutral mixture {Za = 0) left and charged 
mixture {Za 0) right. This allows to inspect changes in temperature behaviours 
of transport coefficients visually (~ T^f^ in Enskog-like theories and ^ T"^/^ at 
presence of long-range Coulomb interparticle interaction). The bulk viscosity n 
does not depend on charge Za and T. In such a way it does not change the nature 
of its temperature behaviours. 

A somewhat other quantitatively pattern appears when Ua np. If weighty 
component of mixture is in greater amount, partial lines in figures become much 
discrete and vice versa. Moreover, at some correlations of densities they can cross 
and the influence of light component becomes predominant. But temperature in- 
creasing diminishes this difference. (More figures see in details on www-server: 
www.icmp.lviv.ua/icmp/preprints/PS/9621Ups.gz). 

At the end of this section one should add a little comment about choosing 
g2^{r) and D, since they appear in transport coefficient expressions. As have been 
obtained earlier in this paper, transport coefficients take an analytical structure 
within the Chapman-Enskog method. But obtaining of any similar relations for 
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§2^ {r) and D, which are perceived as being so adapted for considering system 
closely, is a completely different problem. Thus, in the present paper we have been 
used results done already. The quantity D = l/2r is chosen as a solution to the 
next equation: 



M 



7 ^ 



21- A 



[1 + ra^] 



M 



Pn 



E"-ct< 
, 1 + 



M 

n„a„. 



y- 
^ 1 



2(1 ~ z\) ^ 1 + roo, 

^ ' a— 1 



which has been proposed i n But the expression for g^^ir) has been assembled 
on the basis of results in ||44 |45|, In [Q, the Percus-Yevick equation for 

the radial distribution function in liquid has been generalized to multicomponent 
mixtures, in 92^ {f) has been calculated for ion mixtures at the presence of 
Coulomb and screened Coulomb interactions. A much precise formula for calcula- 
tion of 52 ^ (r) on the basis of the Percus-Yevick approximation has been obtained 
in |^6j, as well as a new closure for the Ornstein-Zernike equation for a system of 
charged particles with specific interparticle interaction potential has been used in 
prf . Theoretical solutions there are compared with Monte Carlo simulations. Fi- 
nally, we conclude that the following structure of g2'^(r) on hard sphere contact is 
the most plausible (in the present version of our work we wish to obtain analytical 
formulae only, choosing actual expressions for §2^: 



92 y^a) 
92 V^afi) 

30L I \ 

92 V^Sia) 



1 + + l^o.<jl{<Tp - a„) \ (1 - A)-' 



'^P92°'{<^a) + (Ja92^{<Jl3) /2c7a(3, 
- (ffa + (Jp) ■ 



Structurally such a configuration is of simple design also. 
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Appendix A. Calculation of A"-coefRcients 



The two mutually exchangeable infinite systems of equations for searching ^"-co- 
efficients could be written in the form 

oo 
oo 

where quantities asr, Ps have an integral structure and are expressed via LJ^ poly- 
nomials. General expressions for these quantities are very cumbersome and their 
forms depend on the structure of collision integrals. But the general form is always 
the same and is completely shown in p^ . In our case 

- J d<d.^,-d. fricmric^.) {^f^L^ (^^) X 



I aa 
92 



2kT 



2kT 



(c?,C?)L3/2 



mg [clY 
2kT 



( mg {err 

ma {c^f 
2kT 



Jab,a) 



'/\3/2(mg (cff 



2kT 



iclOLf/^ ( ^^^^] (c?,c?)L?/^ 1^^^'^^^' 



2kT 



mg (cr) 
2kT 



2kT 



mg {cD 
2kT 



'm, (40' 



2kT 



'm, {cfY 



2kT 



nib 



i4y 



2kT 



mt {4y 



2kT 



13 



2fcT 



2fcT 



2kT 



2n J \ 

1 + mana(7aag2° 



P a, 13 = 1 



Let us consider the case, when indices r and s take two values only: zero and 
unity, which correspond to the first polynomial approximation. In such a case the 
system of equations becomes finite. It is convenient to represent it in the matrix 
form: 



fly Aj 



A, = col [aI Al, Al Aj), k = col (l3i^\ /3f^), 



matrix aij reads 



(ah. a'] (ah, a) 

[ah. a) (ab.a) 

y\o = «2l ail 

(ba.a) (ba.a) 

(6a, a) (ba.a) 

= «4l ail 



: ai2 

■ 0,22 

■ 0.32 

■ 0.42 



(ab,b) 

(at, 6) 

{ba,b) 

Jba,b) 
-*10 



0.13 
023 
033 
043 



(ab,b) 

(at, 6) 
^11 
(6a, 6) 

Jba,b) 



The solution is written using the Kramers' method at once: 



014 
024 
034 
044 



^0 



det cii 



h ai2 ai3 ai4 

62 0.22 023 «24 

63 0.32 fl.33 a34 

64 042 043 044 



Al 



det ai 



an hi ai3 cm 

^21 &2 «23 ^24 

flSl ^3 «33 034 

^41 64 tt43 fl44 



A^ - — 
det< 



an ai2 61 flu 

021 022 ^2 ^24 

031 032 63 034 

O4I fl42 ^4 044 



A\ 



det ( 



On an 013 61 

121 021 023 ^2 

031 031 033 63 

041 041 043 64 



Quantities aij can be calculated and represented via so-called f2-integrals (see Ap 



pendix C). Then one obtains: 



On 



32 m* 



g2°''K6|n,/3)'^''X?;;: 



(1,1) 



(14) 



/32 /m* 

ai2 = 021 = -naTibx — I X 



5aa^(l,l) _ <iQo(l'2) 
2 hs hs 



(1,1) 00^7'^^^^ 



14 



ai4 = 023 = maUbVT^ X 



rub J 



2 hs 



hs 



/32 /to* 

V TT 



52''(^a6|n,/3) 



30 /rrv^y 25 /nr^^^' 
4 l^mfcy 4 Vtoc^ 



aa ^(1,1) _ 



V^a/ V™a/ ^ (TOa + TO6)2 



(2,2) 
hs 



30 /m*\^ 25 /to*" ^ 



4 V TOb / 4 V TOa 



5 -12(1,2) 



(1,3) ^ 2TOaTOb aa^(2,2) 



(TOa + rribY 



+ 



1 (to*)5/2 
024 = 4nan6^ T7TT X 



^afc^a,l) _ r,cib^{l,2) _^ a6^(l,3) _ .,afen(2,2) 



'hs 



hs 



hs 



+ 



^a6^(l,l) _ ^ab^{l,2) ^a6^(l,3) _ 2ab^{2.2) 
^1 1 1 1 

The remaining quantities aij are expressed via previous ones by simple exchange 
of indices a ^ b and vice versa, b ^ a: 





031 


= Oi3 


a- 


b— >a5 


032 


= 041 


= Oi4 


a- 


->6, 6— ^a? 




033 


= On 


a- 


-»6, 6— ^aj 


034 


= 043 


= 012 


a- 


-»6, 6— ^aj 




042 


= 024 


a- 


-»6, 6— ^a? 




O44 


= 022 


a- 


6 — 



The matrix elements bi in the first polynomial approximation are calculated exactly: 

5 / 27r " \ 
bi = -na ( 1 + -^^rnanaal^g2°'{aaa\n,P) ) - 



a=l 
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2 V2 3 , 

\ a=l 



27rmo 

a,/3=l 



5TOa , 



= I a^b, b^a, 

&4 = &2 I a^6, b^a, 



where p = ^ pa = 
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Appendix B. Calculation of S"-coefRcients 

The two mutually exchangeable infinite systems of equations for searching B"-co- 
efficients could be written in the form 



oo 



r=0 



oo 



r=0 



where quantities jsr, Cs have an integral structure and are expressed via poly- 
nomials. General expressions for these quantities are very cumbersome and their 
forms depend on the structure of collision integrals. But general form is always the 
same and is completely shown in feq]. In our case 



d<d^§3'""da;/r^"^(c?)/i"^"'(ci) 



.2kTJ 



2kT 



r5/2 



{cry 

2kT 



,5/2 



-m-a (4'! 
2kT 



r5/2 



5/2 



2kT 

2 



f^i ci - - (ci ) I ] - L 



2kT 



r5/2 



(4) 

2kT 



r5/2 



2kT 



5/2 ( mg (eg*)' 



2kT 



2kT 



C2C2--(C2) / 



^2^2 - - (C2 ) I 



C2C2--(C2) / 



.«) _ / ^„a^„a ab ."(0) , a wKO) , J, s / ™a \2 ,5/2 ( rrig (cf) 



2kT 



at 



5/2 



2kT 



/2 / '"g (c^)' 



2fcT 



r5/2 



{cry 



2kT 



Ci - - (Ci ) / - Lr 



r5/2 



2A;T 



c^c? - - {elf I 



1. 



(ab.b) 



d<dx,i9"'d../,"(°)(c?)/,^('"(c§)^^L^/^ ( 

1 2y Ji V ^ 2;^2fcT)2 I 2fcT 



9?" 



■'" (=#^)(=S'<*' - i ' ) - - i (4)= 7) 



5/2/"X^4c|*)^\/ 1 



2kT 



Cn Cr) — Co 



1^7 



6-1 2 



2kT 



44 --^{4} I 
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2kT 



a — l f 



2^ " 



/3 = 1 " 



I Ci Ci 



47r 
15ma 



J2 ^aUa'^l^aT 



Let us consider the case, when indices r and s take one zero value only, which 
corresponds to the zeroth polynomial approximation. In such a case the system of 
equations becomes finite. It is convenient to represent it in the matrix form: 



Bj — col 



{bs,b^), d, = coi(c^'^), cr), 



matrix Cij reads 



(ab.a) , (a) 

Too + Too = cii 

(ba.a) 

Too = C21 



(ab,b) 

Too 

(ba,b) , (6) 

7oo +7oo 



Cl2 
C22 



The solution is written using the Kramers' method at once: 



r>a 

^0 - det 



1 



dl Ci2 
d2 C22 



^" - det 



1 



cii dl 

C21 d2 



Quantities Cij can be calculated and represented via so-called J7-integrals (see Ap 



pendix C). Then one obtains: 



Cll 



-naTlb- 



3 V TTiria marrib 



gfi'Jab\n,P) 



5aa^(ll)^3m,aa^(2 2) 

2ma 



2ma 



kT 



2) _|_ aa j^(2.2) 



8nani,\/2kTm* 
Cl2 = —, ^ X 



3{ma + nib) 



9f{'7ab\n,/3) 



gafc^Cl,!) _ 3afc„(2,2) 

ns 2 ns 



^ab^ilS) _ 3,,^(2^2) 



The remaining quantities are expressed via previous ones by simple exchange 
of indices a — *■ 6 and vice versa, 5 — > a: 



C21 = Ci2 
C22 = Cll 



a — b — >a 1 
a — >b, b — >a • 
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The matrix elements in the zeroth polynomial approximation are calculated 
exactly: 



a,/3=l 



27r 



3ma 



a=l 



7rn 
+ 



1177rno 



471- 



\ " a=l / 



d2 = di I a^b, b^o- 



19 



Appendix C. !7-integrals 

To calculate J7-integrals analytically, one should find first the expressions for scat- 
tering angles x' and x* ■ In general case both these angles are complicated functions 
of many factors: 

X''* = TT - 26 / dr 1 - ^ 

ro 

here the meaning of ro is searched from the extremum condition 

In view of geometrical reasoning, one can calculate cos (x') exactly (Figure CI): 
cos ("'3x72) = h/a^p. 




Figure CI: Collision of hard spheres. 

In the paper |2^ , cos (x* ) is calculated approximately. In small scattering angles 
approximation one may be thought of as \g'\ ~ \g\ and Ag = Qy. Considering the 
scattering dynamics as is shown in the drawing below (see Figure |C2| ), one may 
express cos (x*) = sin 6' (x* + = 7r/2) via and make expansion in the Teilor 
series. Here F denotes the Coulomb force, which acts on particle 2 from the side of 
particle 1; in considered reference frame particle 1 is fixed while particle 2 impacts 



from infinity. Hence taking into account (3/7) one may obtain 



This allows to write all ]7-integrals in an analytical form, which appear at searching 



^^-coefficients (see Appendix A.), i?"-coefficients (see Appendix B.) and 



20 




Figure C2: Scattering of two charged particles. 



i?" -coefficients: 







a/3 ^(14) 


^3 

2 


V 2fcT 




= 3 




a/3 ^(1:2) 


= 2 


a/3^(l,l)- 


hs 


= 12 


a/3^(l,l)- 
ns 


a/3 ^(1,3) 


= 12 


a/3^(l,l)- 




= 2 


a/3 ^(1,1)' 


a/3 ^(2,2) 


= 4 


a/3^(l,l)- 



where the new quantity D has been introduced, while calculating Q\ . It is 
the cut-off radius of upper limit in the integral with respect to h. This is needed 
to avoid the difficulty of logarithmical divergency in J7;-integrals at long distances. 
The quantity D has dimensionality of distance and meaning of screening radius 
(like a Debye radius). The equation for obtaining D in any explicit form is a 
separate task, which will be considered in a next paper. 
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Figure 1 caption 

The temperature dependencies of some transport coefficients for binary mixtures: 
a), c) and e) for neutral systems at Z\ = 0.1; b), d) and f) for once- ionized gases 
at Z\ = 0.0125. Everywhere T: [K]. a), b) are mutual diffusion coefficients for 
mixtures of Ar-Kr (•),' Ar-Xe (.), Kr-Xe (a), a): 10"^ b): 10-9[mVs]. c), d) 
are thermal diffusion coefScients for Ar-Kr mixture: argon's Ar (+) and krypton's 
Kr (♦) components, _Dt c): 10^'"', d): 10^^[kg/m-s]. e), f) arc shear viscosity 
coefficients for Ar-Kr mixture (•), argon's Ar (.) and krypton's Kr (a) components, 
ri e): 10"^, f): 10-f'[Pa-s]. 
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